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Abstract. We compare the behavior of a wave packet in the presence of a complex and a 
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sheds new light on the possibility to recognize quaternionic deviations from standard 
quantum mechanics. 
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I. INTRODUCTION 

This article represents the third work of the authors on the Schrodinger equation in the presence of 
a quaternionic potential step. In the first paper [1], we have shown that, for such a potential, it is 
possible to calculate an analytic plane wave solution. This represents, to the best of our knowledge, 
the first case in which an analytic solution has been given for a quaternionic quantum mechanical 
potential problem (previous plane wave studies, regarding the barrier [2-4] and well [5] potentials, 
required numerical calculations). The possibility to work with an analytic plane wave solution allowed, 
in our second paper [6] , a detailed discussion of the quaternionic diffusion through the stationary phase 
method. The motivation to write a third paper on this topic is mainly due to the old wish of the authors 
to find qualitative differences between complex and quaternionic formulation of quantum mechanics 
which could be useful in identifying the evidence of quaternionic potentials in the case in which such 
potentials really exist. In this spirit, the quaternionic diffusion, at the moment discussed from a general 
point of view by using the stationary phase method, should be investigated by directly analyzing the 
motion of the incident, reflected and transmitted wave packets in a potential step problem. The 
numerical results of this study can be then interpreted by looking for analytical approximations of 
the incident, reflected and transmitted wave packets in the case of gaussian convolution functions. 
The comparison between the complex and pure quaternionic case for the potential step seems to be 
the best starting point to analyze qualitative differences between quaternionic and complex quantum 
mechanics. Indeed, more complicated potentials can often been seen as successive potential steps. For 
example, a potential barrier can be successful studied as a two-step problem [7, 8] . The investigation 
proposed in this paper could also give a final answer to the old question concerning the possibility 
to fit a pure quaternionic potential by a complex one a consequently to never recognize quaternionic 
deviations from standard quantum mechanics. 



II. PLANE WAVE ANALYSIS 

Before to give the explicit plane wave solution for a pure quaternionic potential step, let us show how, 
by a simple re-phasing of the wave function, we can transform, without lost of generalization, the 
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(j, fc)-part of the quaternionic potential 



iV.+jV.+kV, , 



into a pure j-part, i.e. 

This will greatly simplify our presentation. Let us begin by considering the quaternionic Schrodingcr 
equation [9] 

ft 2 

i 1 ^-^xx{x,t)-(iV 1 +jV 2 + kV 3 ) $(x,t) = h$ t (x,t) . (1) 
Am 

This partial differential equation can be reduced, by using the well-know separation of variables 

t) = ip(x) exp[— i Et/h] , 
to the following second order ordinary differential equation [10] 



i — <p"(x) - (i K + jV 2 + kV 3 ) <p(x) = - <p(x)Ei 
2m 



(2) 



By multiplying (from the left) the previous equation by an unitary complex number e la and re- writing 
the pure quaternionic part of the potential as 



jV 2 + kV 3 =j^V 2 2 +V 3 2 e 



where 9 = arctan [V^/V^], we find 

h 



2m 



^(x)-^V 1+ j^/v:+V 3 2 e-^ 



(p(x) 



■e ia ip(x)Ei 



By observing that e la j = j e la , we can rewrite the previous equation as follows 



-i(a+9) 



ip{x) 



f{x)\ E ! 



(3) 



(4) 



(5) 



Consequently, the choice a — —9/2 conduces to 
where 

ip(x) = e~ l z (p{x) . 

It is important to observe that constant phases paly no role in the stationary phase method and, thus, 
such a re-phasing has no physical effect on the motion of the incident, reflected and transmitted wave 
packets. For the diffusion problem, £ > y o , wc shall compare the complex case, V = V t (V 2 = V 3 = 0), 

with the pure quaternionic case, V = \/V 2 2 + V 2 (Vi = 0). 



II.A COMPLEX POTENTIAL STEP 

By setting V 1 = V and V 2 = V 3 = in Eq.(5), we find the standard Schrodinger equation, 

ft 2 

i — if>"(x) - i V ip c {x) = - ip c (x) E i , 
2m 

whose analytic plane wave solution is given by [11] 

[I] free region (x < 0) : ip C)I (e,x) = exp[ i e x } + r c (e) exp[ - i e x ] 



[II] potential region (x > 0) : ip c>II (e,x) = t c (e) exp[iax] , 



where 
and 



e = V2m E/H , a = v / 2m (E - V„)/h , 
r c (e) = (e - a)/(e + a) , t c (e) = 2 e/(e + a) . 



(0) 



(7) 
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II. B PURE QUATERNIONIC POTENTIAL STEP 



(9) 



By setting in Eq.(5) V 1 = and \JV 2 + V 2 =V , we obtain 

ft 2 

i — ^{x)-jV Q %{x) = -^ q {x)Ei. (8) 

The analytic plane wave solution for the pure quaternionic potential step reads 

[I] free REGION : ip qjI (e,x) = exp[ i e x ] + r q (e) exp[ — i e x ] + j r q (e) exp[ e x ] , 

[II] potential region : i/jg itI (e,x) = (1 + j w) t q (e) exp[i p x] + (w + j) t q (e) exp[- px] 
where 

e = V2mE/h , p = ^2m^E 2 -^/h , w = - i V / (e + \J E 2 - V*\ , 

and 

r q (e) = (e - p) cxp [i arctan (e/p) ] / ^e 2 + p 2 , 
f g (e) = (l + i)ew/{e + p) , 
= e /P ; 

t ? (e) = Ve 2 + P 2 e w exp [ - i arctan (p/e) ] / [p (e + p)] . 

For a detailed derivation of the plane wave solution for a quaternionic potential step, we refer the 
reader to the paper cited in rcf. [1]. 



III. WAVE PACKET ANALYSIS 

Until now, we have been concerned only with plane waves. In this section, we are going to study the 
time evolution of quaternionic wave packets and deducing from them several important results. The 
principle of superposition guarantees that every real linear combination of the plane waves in region 
I and region II, 

® ItII (x,t)= Jdeg(e)^ ItII (e,x)eiq ? [-iEt/h} [E = {he) 2 / 2 m] , 

will satisfy the Schrodinger equation in the presence of a potential step. The reason why the use of 
the wave packet formalism is very interesting lies in the fact that the calculations, in the case of a 
gaussian convolution function 

g(e) = cxp [a 2 (e~e 2 )/4}/2^ , 

can be analytically approximated. This will allow to check and interpret our numerical results. In 
the next subsections, we first discuss the wave packets motion in the case of a standard (complex) 
potential step and, then, we analyze the special case of a pure quaternionic potential. 



III.A COMPLEX CASE 

The wave packets in region I and region II are given by 
Vt CJ {x,t) 



where 



/■+ oo 

/ d(a e) g(e) { exp[i e x ] + r c {e) exp[ — i e x ] } exp[ — i e 2 htj 2m] , 

J a e m i n 
r+ oo 

> u(x, t) = / d(a e) g(e) { t c (e) exp[i a x ] } exp[ — i e 2 htj 2m] , 

J a e m i n 



a£ min = a\/2m V / h 
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We find three wave packets: incident, reflected and transmitted, 



fi c ,taoOM) 



+ oo 



d(ae) (?(e) cxp[z (ei-e 2 ht/ 2m) ] 



+ OO 



d(a e) g(e) r c (e) exp[ — i {ex -\- e 2 ht/ 2m) ] 



+ OO 



d(ae) g(e) i c (e) exp[z (c 



<ht/2m)] 



The choice of a gaussian convolution function g(e) peaked in e (e > e min ) and whose value is 
practically zero near to e mi „, i.e. g(e min ) ~ 0, allows to legitimately approximate the incident wave 
packet as follows 



fi C)lnc (:r,i) 



+ OO 



d(a e) g(e) exp[i (ex — e 2 fit/ 2m) ] 



cxp[i (e x — a 2 e 2 r/2) — i arctan(2r)/2] 



exp [— (x/a — ae r) 2 / (1 + 2 i r)] 



(1 +4r 



vl/4 



(10) 



where r = ht/ma 2 . In the case in which the variation of r c (e) can be neglected compared to that of 
c/(e), the reflected wave has the same form as the incident wave packet [11], 



fi c ,rcf(z,t) ~ r c (e ) / d(ae) g(e) exp[- i (ex + e 2 ht/ 2m) ] 

./ a e m i n 

= r c (e )Q. Ctino (-x,t) . 



(11) 



For the transmitted wave packet, by neglecting the variation of t c (e) compared to that of g(e) and by 
approximating a by its first order Taylor expansion, a + (e — e ) e /er , we find 



n c ,tra(a;,t) ~ ic(eo) exp[z (cr 2 - e 2 ) x/cr ] / d(ae) p(e) exp[i (ee x/a - e 2 ht/ 2m) ] 

= t c (e ) exp[z (cr 2 - e 2 )x/cr ] O c , illc (e x/a ,t) . (12) 

The reflection probability (the ratio between the probabilities of finding the particle in the reflected 
packet, at a positive time t , and in the incident packet, at time — 1 ) is given by 



dx | Q Cjrcf (x,i ) | 



dx | f2 C;inc (x, - 1 ) 

OO 

•0 



|r c (e )| / dx\n Ctinc (-x,t 
\r c (e )\ 2 



dx I n Citoo (ar, - < ) 



(13) 



Similarly, the transmission probability is given by 



drc | n c ,t™(a;,*o) I / / cLc|n C)inc (a;,-to) | 

J — oo 

2 i r° 

I Citac (e a;/<7 0! t ) | / / dx | ^ c , inc (a;, - t ) | 

J — OO 



2 / ,+ °° 

|*c(e )| / dx 
Jo 



K( £ o)| 2 - 



c o J 



dx I fi cloc (x,i ) 



dx | r2 Cjinc (x, — i 



- \tc(e )f 

Co 



(14) 



It is easy to verify that 



,t 



4e 2 



e (e + a) 2 



= 1 . 
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III.B PURE QUATERNIONIC CASE 

The wave packets in region I and region II are given by 

Q q>I (x, t) = d(a e) g(e) { exp[i ex] + r q (e) exp[ — iex] + j r q (e) exp[ ex] } exp[ — ie 2 htj 2m] , 

J a e m i n 
p+ oo 

Qq,n{x, t) = d(a e) g(e) { t q (e) exp[i px] +wt 9 (e) exp[ - pi] } cxp[ — ie 2 htj 2m ] + 

J a e m i n 
/> oo 

j / d(a e) g(e) { w t q (e) exp[i px ] + t q (e) exp[ — px ] } cxp[ — i e 2 htj 2m ] . 

J a € m i n 

For sufficiently large negative t, only the (complex) incident wave packet 

r+ oo 

ft q ,inc(x, t >C 0) = / d(ae) g(e) expfi (ex - e 2 ht/ 2m) ] — Q Ciinc (x, t) 



exists. For sufficiently large positive t, the evanescent wave packets disappear and we only find one 
(complex) reflected wave packet, 

r+ oo 

^ q ,vat(x,t ^> 0) — / d(a e) g(e) r q (e) exp[— i (e x + e 2 ht/ 2m) ] w r q (e ) O c ,mc(— x, t) , 

J a e m i n 

and one (quaternionic) transmitted wave packet, 

r+ oo 

n q ^(x,t > 0) = / d(ae) g(e) (1 + jw)t q (e) cxp[i(px- e 2 ht/2m)] 

J a € m in 

~ (1 + j w ) t q (e„) exp[i(p* -<%)x/p 3 ] n c , inc (e 3 x/p 3 , t) . 
A simple calculation shows that 



P q , tBf = dx\Q q , tBf (x,t )\ 2 j / dx\n q ^ c (x,~t )\ 2 » |r g (e )| 2 (15) 

J — oo J — oo 

and 

/+ OO /• 3 . . 

dx|l\ tra (x,t )| 2 / / da;|fi g , lnc (x,-t )| 2 w J (l + W'j |i g (e )| 2 . (16) 



Finally, 



Pq,r<si + Pi J,tra 



- />) , P 3 ^ x , £ 2 -P 2 \ e 2 



+ p* e \ e + P 



1 



IV. NUMERICAL RESULTS AND INTERPRETATION 

The results of our numerical study are plotted in Fig. 1, where the probability densities |f2 c (x,t)| 
(complex case) and |f2 9 (x,i)| (pure quaternionic case) are drawn as a function of x/a for different 
values of r = ht/ma 2 . The incident wave packets are centered in E = 2 V (diffusion) and the potential 
is given by a y/2m V„ = 10 2 h. The data show an interesting phenomenon. For the same potential and 
the same incident energy, the reflected waves have different amplitudes in the case of a complex or 
a pure quaternionic potential step, and the transmitted waves move with different velocities. The 
quaternionic transmitted wave packet moves faster than the complex one. Let us try to understand 
these results by using the analytic approximations given in section III for the incident, reflected and 
transmitted wave packets and by using the fact that for small times, t < 1, no spreading effect is 
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present. This allows to simplify the expressions given in section III for the complex and quaternionic 
wave packets. The amplitude of the incident wave is then given by 



\tt q ,inc(x, t)\ w |O c , loc (x, t)\ a cxp [- (x/ 



a — a e„T) 



For the reflected waves, we have 



|f2 c ,refOM)| ~ kc(e )| exp [- (x/a+ ae r) 2 
\Qg ttet (x,t)\ » |r g (e )| cxp [- (x/a + ae r) 2 



This implies 



ic;>,t)i7ic:i(x,t)i 

|^ rof (x, £)| 2 / |f\ rof (x,t)| 2 



|r c (e )| 
|r- 9 (e )| 2 

k c (e )| 2 / |r ? ( £o )| 2 



2.94 x 10" 2 
2.58 x 10" 3 
11.4 , 



which confirms the numerical results given in Fig.l (negative x-axis). We can also see that the reflected 
waves move with the same velocity, 



(17) 



For the transmitted waves, by observing that 

|fi C)tr a(a;,t)| » |* c (e )| cxp [- (e /er ) 2 (x/a- aa a r) 2 ] , 
|fi 9 ,tr»(aJ,t)| ~ 1(1 + jw a )t q {e a )\ exp [- (e /p ) B {x/a - aplr/e 2 ) 2 ] , 

we find that the centers of the complex and quaternionic transmitted wave packets move, respectively, 
with velocities 

w c ,tra = K 0"o / m and «q,tra = ^-Po/ m e 2 . 

Consequently, 



(18) 



- 1 



Eo 

v 



a\/2m V Ht 



ma- 



V 



1/4 



and by using 



we find 



E = 2V 



a\J2m V 
H 

a 



= 10 

: lOV 
! 1.14 . 



ht 



ma 



This agrees with our numerical calculations as shown by the motion of the complex and quaternionic 
transmitted wave packets plotted in Fig.l (positive x-axis). For example, at time t = 0.15 the max- 
imum of the complex transmitted wave is found at x = 15 a, whereas the center of the quaternionic 
wave packet (which moves faster) reaches at the same time the point x ~ 17.1 a. 



V. CONCLUSIONS 

In the last years, the Schrodingcr equation in the presence of quaternionic potentials with constant 
and space-dependent phase has been a matter of study and discussion in the literature [4-6]. Some 
properties of this class of potentials are detailed discussed in the Adler book [9] which represents a 
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milestone in such a field. If quaternionic quantum mechanics represents a possible way to describe the 
nature, then it becomes relevant to examine how the predictions of standard theories may be affected 
by changing from complex to quaternionic potentials. The first theoretical analysis of quaternionic 
potential barriers was developped in [2, 3] and showed that in contrast with the standard complex 
case where the left and right transmission probability are equal in magnitude and phase, in the 
quaternionic quantum mechanics only the magnitude are equal. So, the measurement of a phase shift 
should be an indicator of quaternionic effects and space-dependent phase. Nevertheless, as remarked 
in [9], experiments to detect a phase shift are equivalent to detecting time-reversal violation and 
consequently cannot be seen in neutron-optical experiments [12-14]. A more complete phenomenology 
of the quaternionic potential barrier is given in [4] . With respect to previous works regarding potential 
barrier diffusion, in this paper, we have introduced the quaternionic wave packet formalism. This allow 
to study new qualitative differences between standard quantum mechanics and theoretical solutions 
obtained by solving the Schrodinger equation in the presence of a quaternionic step. For a detailed 
discussion of experimental proposals on quaternionic potentials and quantum mechanical systems in 
which a CP violation is interpreted by potentials with a space-dependent phase, we refer the reader 
to references [1,4]. 

In this paper, we have seen that for incoming particles with an energy spectrum centered in 
E a = 2 V the complex and quaternionic transmitted wave packets move with different velocities, 



u q,ti. 



V3V3 



Due to the fact that the complex and quaternionic incident wave packets coincide for large negative 
times (incoming particles), we could try to fit the motion of the quaternionic wave packet by using 
a different complex potential step. For example by choosing a new complex potential step, W , such 
that 

-,3/4 



^-1 = 

W 



Eo 
V 



=M -i 



we find 



</.t. 



>(*) 



^(t). In the particular case investigated in this paper, E = 2 V , this implies 



3V3 



W 4 

In this scenario, we have the same incident wave and complex or quaternionic transmitted wave 
packets which move with the same velocity (note that the reflected waves, for a same incoming 
energy spectrum, always move with the same velocity because propagate in the free potential region). 
Once we have guaranteed same velocities in the potential region, we have to check the probability 
of transmission (or equivalently the probability of reflection). A simple calculation (see section III) 
shows that 



^K-VEo- Wo 



rWo-^JEo-Wo 



E n - v: 



1/4 



E a - V E„ - V„ 



2.01 x 10" 



2.58 x 10" 



This gives an interesting answer to the old question concerning the possibility to fit a quaternionic 
potential by a complex one and consequently to never recognize deviations from standard quantum 
mechanics. The study presented for a very simple potential, i.e. the potential step, shows that incoming 
particles with a given energy spectrum behave differently if they are diffused by a complex or a pure 
quaternionic potential step. This result, obtained by numerical calculations and interpreted by our 
gaussian analytical approximations, cannot be seen in the plane wave analysis. In such a limit, we only 
have reflection and transmission probabilities, and the quaternionic case can always be interpreted 
in terms of a complex potential step which appropriately fits the quaternionic probabilities. This 
paper stimulates further investigations, for example it should be interesting to study, within the wave 
packet formalism, diffusion and tunnelling by complex and pure quaternionic potential barrier which, 
surely, represents a more realistic potential to be tested by an experimental analysis. In the case 
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of the potential barrier, it is of great interest to examine the tunnelling phenomena, in particular, 
in a forthcoming paper we aim to analyze in detail the Hartman effect [15, 16] in the presence of 
a pure quaternioic potential. The possibility to treat the potential barrier problem as a two-step 
problem [7, 8] gives a more solid interest to this paper, what it could appear an academic quantum 
mechanical exercise play a fundamental role in the understanding the behavior of wave packets in 
the presence of a more complicate potential. The analysis done in this paper explicitly shows that a 
qualitative difference exists between complex and quaternionic quantum mechanics. 
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Fig. 1. Diffusion [E = 2 Vo) of a wave packet at a complex and pure quaternionic potential step. The complex 

(Vo = Vi) and quaternionic (Vo = \J V% + V% ) potentials act for positive x. For negative time (r < 0) only 
the incident packet is present and it moves towards the step. After a certain time, we find four packets. The 
reflected packets, \Q Ctinc (x,t)\ 2 and |f2 9lill c(a;,t)| 2 , are returning to the left moving with the same velocity. 
The transmitted packets, |f2 c ,tra(£,t)| 2 and \Q. q ,ti&{x,t)\ 2 , propagate towards the right and move with different 
velocities. 



